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We solve the 3-loop ^-derivable approximation to the ther- 
modynamics of the massless field theory by reducing it to a 
1-parameter variational problem. The thermodynamic potential is 
expanded in powers of and m/T, where g is the coupling con- 
stant, m is a variational mass parameter, and T is the temperature. 
There are ultraviolet divergences beginning at 6th order in g that 
cannot be removed by renormalization. However the finite thermo- 
dynamic potential obtained by truncating after terms of 5th order 
in g and m/T defines a stable approximation to the thermodynamic 
functions. 

The thermodynamic functions for massless relativis- 
tic field theories at high temperature T can be calculated 
as weak-coupling expansions in the coupling constant g. 
They have been calculated explicitly through order 
for the massless field theory for QED and 
for QCD dlJ^. Unless the couphng constant is tiny, the 
weak-coupling expansions are poorly convergent and sen- 
sitive to the renormalization scale. This makes the weak- 
coupling expansion essentially useless as a quantitative 
tool: it seems to be reliable only when the coupling con- 
stant is so small that the corrections to ideal gas behavior 
are negligibly small. The physical origin of the instability 
seems to be effects associated with screening and quasi- 
particles. 

A possible solution to this instability problem is to 
reorganize the weak-coupling expansion within a varia- 
tional framework. A variational approximation can be 
defined by a thermodynamic potential £7 that depends on 
a set of variational parameters m^. The free energy and 
other thermodynamic functions are given by the values 
of U, and its derivatives at the variational point where 
dfl/dnii = 0. A variational approximation is system- 
atically improvable if there is a sequence of successive 
approximations to Q that reproduce the weak-coupling 
expansions of the thermodynamic functions to succes- 
sively higher orders in g. An example of a systematically 
improvable variational approximation is screened pertur- 
bation theory, which involves a single variational mass 
parameter |^]. 

A variational approach can be useful only if the cor- 
rect physics can be captured by appropriate choices of 
the variational parameters. Information about screening 
and quasiparticle effects is contained within the exact 
propagator of the field theory. The possibility that these 
effects are responsible for the instability of the weak- 
coupling expansion suggests the use of the propagator 
as a variational function. Such a variational formulation 



was constructed for nonrelativistic fermions by Luttinger 
and Ward 1^ and by Baym |^ long ago. In the case of 
a relativistic scalar field theory, the propagator has the 
form [P"^ + Il{P)]~^ , where n(P) is the self-energy which 
depends on the momentum P. The thermodynamic po- 
tential has the form 



^^o[n] = i^ [iog(p2 + n) 
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where the interaction functional $[11] can be expressed 
as a sum of 2-particle-irreducible diagrams. It is con- 
structed so that the solution to the variational equation 
SHq/SH = is the exact self-energy, and the value of flo 
at the variational point is the exact free energy. We can 
obtain a systematically improvable variational approxi- 
mation by truncating $ at n'th order in the loop expan- 
sion, where n = 2,3, .... We refer to such an approxima- 
tion as the n-loop ^-derivable approximation. The 2-loop 
<I>-derivable approximation for QCD has recently been 
used as the basis for quasiparticle models of the ther- 
modynamics of the quark-gluon plasma pC| , pT| . Since 
(^-derivable approximations guarantee consistency with 
conservation laws, they may be particularly useful for 
nonequilibrium problems . 

While the ^-derivable approximation is easily for- 
mulated, it is not so easy to solve. The variational equa- 
tion is a nontrivial integral equation that, to the best 
of our knowlege, has never been solved for a relativistic 
field theory except in trivial cases where the self-energy 
is independent of the momentum. The main problem is 
that the thermodynamic potential has severe ultraviolet 
divergences that vanish at the variational point only if 
<I> is calculated to all orders. They do not vanish away 
from the variational point, and they do not vanish at the 
variational point if the loop expansion for $ is truncated. 

In this Letter, we solve the 3-loop ^-derivable ap- 
proximation for a massless scalar field theory with a 
(p* interaction by reducing it to a 1-parameter varia- 
tional problem involving a mass parameter m. The re- 
sulting thermodynamic functions have ultraviolet diver- 
gences beginning at 6th order in the coupling constant g. 
However the finite thermodynamic potential obtained by 
truncating after terms of 5th order in g and m/T defines 
a stable approximation to the thermodynamic functions. 

The lagrangian for the massless scalar field theory 
with a (j)'^ interaction with bare coupling constant ao ~ 
5oV(4^)' is 
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FIG. 1. Free energy divided by that of the ideal gas for 
the weak-couphng expansion (dashed lines) and the truncated 
"I>-derivable approximation (solid lines). 

We define the renormalized coupling constant a = 
(7^/(47r)^ using dimensional regularization in 4 — 2e di- 
mensions and the modified minimal subtraction (MS) 
prescription with renormalization scale /i: 



ao/i '^'^ = a 
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The two-loop beta function for the running coupling con- 
stant a(/x) is /3(a) = 3a^ — ^ot'- The thermodynamic 
functions for this theory are known to order ■ The 

weak-coupling expansion for the free energy density is 



= 1 + 15 
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where a — a(/i), L = log(/i/47rr), and J-idcai is the pres- 
sure of the ideal gas of a free massless boson: J^dcai — 
-(7rV90)r4. 

The instability of the weak-coupling expansion is il- 
lustrated in Fig. 1, which shows the free energy divided 
by that of the ideal gas as a function of g{2TrT). The 
dashed lines are the predictions of the weak coupling ex- 
pansion (^ with ^ = 2'itT truncated after orders 5" for 
n = 2, 3, 4, 5. The dashed line for n = 2 is hidden under 
the solid line labelled . The successive approximations 
show no sign of converging. 

To illustrate our method in the simplest possible 
context, we first consider the 2-loop ^-derivable approx- 
imation. The self-energy is independent of P, so we de- 



note it by n 
is 



The interaction functional in rtl) 
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The measure for the dimensionally regularized sum- 
integrals in (|l|) and (||) includes a factor of {e^ ^j? / AttY ^ 
where 7 is Euler's constant. We insert a factor of 
ix~^^ into the sum-integral in (|l|), so fig is independent 
of the renormalization scale. The variational equation 
d^a/dm ~ reduces to the simple gap equation 
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The sum-integrals in (|^) and (^) contain ultraviolet di- 
vergences. Some of the divergences vanish at the vari- 
ational point. They can be removed by adding a term 
proportional to the square of the gap equation: 



no - (m^ - G)V(327r2ao) 



(7) 



where G is the expression on the right side of (|^). The 
additional term gives no contribution to those thermo- 
dynamic functions that are determined by the value of f2 
and its first derivatives at the variational point. The 
remaining divergences can be absorbed into a renor- 
malized coupling constant q;(^) defined by ag/i^^' = 
a{l~a/e)^^. Expanding in powers of m/T, the resulting 
thermodynamic potential f2 is 



^^/^idoai = 1 + 15 [3m^/a - + 4j^3 
+3(L + 7)m^ + ., 



(8) 



where m = m/(27rT). The coupling constant a runs 
with the beta function /9(a) = a^, whose first coefficient 
is too small by a factor of 3 compared to that of the true 
coupling constant. Thus a can be identified with the 
true coupling constant a only at a single scale /xq. If one 
expresses the thermodynamic potential in terms of the 
true coupling constant a defined by (||), the ultraviolet 
divergences cannot be eliminated. The closest one can 
come to defining a finite thermodynamic potential is to 
truncate (0) after the term of order rir^: 



f^/-?^idcai = 1 + 15 [3m^/a — 



4m'' 



(9) 



The gap equation obtained by varying (^) with respect 
to m with a fixed is 



= a (i — m) 



(10) 



Solving this quadratic equation and substituting the so- 
lution into (^, the free energy is 

^/^idcai = 1 - 5a/4 [1 + 6a + 6a^ 



-4 (2 + 3a) vW6 + aV4 , (11) 

where a = a(^). When expanded in powers of g, this 
agrees with the weak-coupling expansion (0) through 
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order . It depends on the renormalization scale /i 
through the dependence of a(/Lt) on /x. However, it is 
much less sensitive to variations in /j, than the truncated 
weak-coupling expansion. The reason is that the trun- 
cated weak-coupling expansion grows like a power of a 
in the strong-coupling limit, while the expression ( p^ ) 
approaches the limiting value |i. 

We now proceed to consider the three-loop $- 
derivable approximation. The interaction term in is 
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where Q4 = — (Q1 + Q2 + Q3) hi the definition of the bas- 
ketball sum-integral. The variational equation obtained 
by varying with respect to n(P) is 



n(p) 



1„2 -2£c> , 
2.9oM ^tad 



1 4 — 4ec^ 



(13) 



^sun 



^p)=i n 

^QiQ2 



Qf + n(Q,)' 



where Q3 = ~(^' + Qi + Q2) hi the definition of the 
sunset sum-integral. The variational equation (13) is a 
nontrivial integral equation for n(P) whose solution is 
complicated by the presence of severe ultraviolet diver- 
gences in the sum-integrals. 

Our strategy is to introduce a variational mass pa- 
rameter m that is of order gT in the weak-coupling limit 
and calculate the sum-integrals as double expansions in 
g and m/T. The variational equation ( [l^ ) is then solved 
for n(P) as a function of P and m. Inserting the solu- 
tion into (^ and ( [l^ ) and expanding in powers of g and 
m/T^ the thermodynamic potential f2 reduces to an al- 
gebraic function of m. The final step is to minimize 
with respect to the variational parameter to. 

We define the mass parameter to implicitly by the 
equation 



^2 _ lf,2,,-2eo. 



1 4 — 4ef^ 



(0,p) . . (14) 



The solution to this gap equation for to can be inter- 
preted as the screening mass. There are two important 
momentum scales: the hard scale 27rT and the soft scale 
TO. The hard region for the momentum P — {2TinT^ p) 
includes n 7^ for all p and also n = with p of order T. 
The soft region is n = and p of order to. We will solve 
the variational equation in the two momentum regions 
separately. We expand n(P) in powers of g^ and m/T. 
For hard momentum P, the expansion has the form 



n(P) = to2 + [n4,o(P) + n4.i(P) 



-^■■-^^ [n8,_2(P) 



(15) 



where n„^fc(P) is of order T'^{m/TY when P is of order 
T. For soft momentum P = (0,p), the function (j{p) in 



the self-energy 11(0, p) 



<j{p) has the form 



o-(p) = .9o/" K -2(p) + 0-4,0 (p) + • . 



(16) 



where an,kip) is of order m?{m/T)^ when p is of order 
TO. We insert the expansions (|l|) and ([l|) into the vari- 
ational equation ( p^ ) and expand in powers of g^ and 
to/T. Matching the coefHcients at each order in g"^ and 
m/T, we find a recursive structure in which the functions 
^n,k{P) and iJn,k{p) are either completely determined or 
expressed in terms of lower order functions. For exam- 
ple, the solutions to the first few self-energy functions at 
hard momentum P are 



n4,o(P)--i^^^^^^^ + ^T^/. 
n4,i(P) = -\Th'}j. 



,(ito), (17) 



Q2(F + g)2 (g2) 



where = —{P + Qi + Q2)- The solutions to the first 
few self-energy functions at soft momentum (0, p) are 



0-4,-2 (p) = -^P^ [Isun{p) - Isun{im) 





(18) 
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The sunset integral appearing in ( [l7| ) and (|lj) is 



Isnn{p) 



^ 1 
qiq2 '?^ + ^ 



where qa = -(p + qi + q2). 

Having solved the gap equation, we can reduce the 
bare thermodynamic potential to a function of the 
single variational parameter to by inserting the expan- 
sions ( p^ ) and (|l^) into (|^) and (|l^) and expanding in 
powers of and m/T. The function f^o contains ultra- 
violet divergences in the form of poles in e. Some of them 
are eliminated when is evaluated at the solution of the 
gap equation. We can cancel them through 5tli order in 
g and m/T by adding a term proportional to (to^ — G)^ , 
where G is the expression on the right side of the gap 
equation (p4|): 
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(19) 



We have introduced an arbitrary momentum scale A in 
order that all terms have the correct dimension 4 — 2e 
when e 0. Other ultraviolet divergences in flo can be 
removed by using (j^) to renormalize the coupling con- 
stant. There are still other ultraviolet divergences at 6th 
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and higher orders in g and m/T that cannot be canceled 
and represent unavoidable ambiguities in the ^-derivable 
approximation. The most severe divergences at 6th or- 
der are double poles in e proportional to am^, a^vn? and 
. Expanding ( |l9| ) in powers of g and m/T and then 
truncating after terms of 5th order, we obtain a finite 
thermodynamic potential: 

(ri/J^idoai - 1) = im^/a 

+ [-m^ + Arh^ + 3 (3L - 2^ + 7) m^] 
+a [2 {I -2 log m - 2.757642) 

-12(^ + 7)77??] 
W [- i - 4 log m - 9.873296) 

+2(^-21og2 + 7)m] , (20) 

where a — a(^), rh — m/(27rT), L — log(/i/47rT), and 
^ = log(A/47rT). 

The gap equation obtained by varying the finite ther- 
modynamic potential (j2^) with respect to with ji and 
A fixed is 



a [i - m - (3i - 2£ + 7) fh^] 
W [-i(^-21ogTO- 3.75764) 



+3 {I + 7) rh] 
[-j^-i(^-21og2 + 7)m] /w? 



(21) 



If we solve this gap equation itcratively in powers of g 
and insert the solution into (20), we recover the weak- 
coupling expansion (jj) for the free energy. Solving the 
gap equation numerically, we find a solution only below 
a critical value of g that depends on /i and A. For /i = 
A = 27rT, the critical value is 5 = 2.61. For larger values 
of g, the thermodynamic potential (|2^) has a run- away 
minimum at m = 0. A similar behavior was observed 
in screened perturbation theory at 3 loops When 
the screening mass is used as the mass parameter, the 
solution to the gap equation terminates at nearly the 
same critical value of g. 

By truncating the thermodynamic potential ( |20| ) 
and the gap equation ( pT] ) after terms of n'th order in 
g and m/T, where n = 2,3,4,5, we obtain a series of 
successive approximations to the free energy. The sta- 
bility of these successive approximations is exhibited in 
Fig. 1, which shows the free energy divided by that of 
the ideal g function of g{2TTT) for ^ = A = 2-kT. 

For the g"^ and g^ truncations, the solutions to the gap 
equation terminate at critical values of g. In contrast to 
the weak-coupling expansion, the predictions of the trun- 
cated ^-derivable approximation seem to be converging 
for values of g below these critical values. The truncated 
^-derivable approximations depend on the renormaliza- 
tion scale ^ and also for n = 4, 5 on the momentum scale 
A introduced in (p^. However the changes in the pre- 
dictions from varying the arbitrary scales is significantly 
smaller than in the weak-coupling expansion. 



We have solved the 3-loop ^-derivable approxima- 
tion for the thermodynamics of the (/)^ field theory by 
reducing it to a problem with a single variational param- 
eter m. We constructed a finite thermodynamic poten- 
tial by adding a term proportional to the square of the 
gap equation and truncating after terms of 5th order in 
g and m/T. This thermodynamic potential has a min- 
imum as a function of m only for g below some critical 
value. Below this critical value of g, the successive ap- 
proximations obtained by truncating the thermodynamic 
potential at orders n = 2,3,4,5 give predictions for the 
pressure that are stable with respect to both the order 
of truncation and variations in ^ and A. The resulting 
predictions for the pressure are numerically close to those 
of screened perturbation theory |13|. The advantage of 
our <I>-derivable approximation is that it represents an 
infinite-parameter variational approximation. Our solu- 
tion to the ^-derivable approximation can in principle be 
extended to higher orders in the loop expansion. Thus it 
provides a systematically improvable approximation to 
the thermodynamic functions that seems to have very 
good convergence properties. 
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